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Abstract 

A d-dimensional finite quantum model system confined to a general hypercubical 
geometry with linear spatial size L and "temporal size" 1/T ( T - temperature of 
the system) is considered in the spherical approximation under periodic boundary 
conditions. Because of its close relation with the system of quantum rotors it rep- 
resents an effective model for studying the low-temperature behaviour of quantum 
Heisenberg antiferromagnets. Close to the zero-temperature quantum critical point 
the ideas of finite-size scaling are used for studying the critical behaviour of the 
model. For a film geometry in different space dimensions \a < d < |cr , where 
< a < 2 controls the long-ranginess of the interactions, an analysis of the free 
energy and the Casimir forces is given. 



1 Introduction and brief overview of the problem 



According to the present understanding the Casimir effect is a phenomenon 
common to all systems characterized by fluctuating quantities on which some 
external boundary conditions are imposed [1,2]. Casimir forces arise from the 
influence of one portion of the system on the fluctuations that occur in a 
portion some distance away. 

The Casimir force in statistical-mechanical systems is usually characterized 
by the excess free energy due to the finite-size contributions to the free energy 
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of the system. The simplest and the most often case of interest is that one of 
a film geometry. 



A simple 0(n) symmetric n > 1 system with a geometry L x oo 2 (and under 
given boundary conditions r imposed across the direction L) is a standard 
statistical mechanical model for describing a fluid, or a magnet, confined be- 
tween two parallel plates of infinite area. One important quantity which arises 
naturally in the thermodynamics of these confined systems is 

(1) FJ... ) {T,L) = - 9 J^R, 

where f^ x (T, L) is the excess free energy 

(2) f?(T,L) = f T (T,L)-Lf bulk (T) 

Here f T (T,L) is the full free energy per unit area (and per k B T) of such a 
system and fbuik(T) is the bulk free energy density. 

In the case of a fluid (then one actually has to consider the excess grand 
potential per unit area and the derivative is performed at constant chemi- 
cal potential \i and temperature T) F(...) is termed the solvation force [3,4] 
F (solvation) i (sometimes called also the disjoining pressure) whereas in the case 
of a magnet one speaks instead about the Casimir force [1,5] F(casimir) (where 
the derivative is performed at constant temperature and magnetic field H). 
In the remainder we will use only the term Casimir force. 

According to the definition given by Eq. (1) the Casimir force is a generalized 
force conjugated to the distance between the surfaces bounding the system 
with the general property F(Casimir)(T, L) — > 0, when L — > oo. We will be 
interested in the behaviour of F( Ca simir) when L >> 1, which is a condition 
for the applicability of the finite-size scaling theory. In general the sign of 
the Casimir force is of particular interest. It is believed that if the boundary 
conditions r are the same at the both surfaces, F^casimir) will be negative (see, 
e.g., [4,6]; strictly speaking, for an Ising-like system this should happen above 
the wetting transition temperature T w [4,6,7]). For boundary conditions that 
are not identical at both confining the system surface planes the Casimir force 
is expected to be positive [4,6]. For a fluid confined between walls this implies 
that the net force between the plates for large separations will be attractive in 
the former and repulsive in the last case. 

In general, the full free energy of a rf-dimensional critical system in the form 
of a film with thickness L, area A, and boundary conditions a and b on the 
two surfaces, at the bulk critical point T c , has the asymptotic form 

(3) U(T C , L) - Lf bulk (T c ) + f: urface (T c ) + f b surface (T c ) + L~^ d ~^A ajb + ■■■ 
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in the limits A — > oo, L >> 1. Here f surface. ls the surface free energy contribu- 
tion and A a 6 is the amplitude of the Casimir interaction. The L dependence 
of the Casimir term (the last one in Eq. (3)) follows from the scale invariance 
of the free energy and has been derived by Fisher and de Gennes [8]. The 
amplitude A a fe is universal, depending on the bulk universality class and the 
universality classes of the boundary conditions [1,5]. 

Equation (3) is valid for both fluid and magnetic systems at criticality. Promi- 
nent examples are, e.g., one-component fluid at the liquid-vapour critical 
point, the binary fluid at the consolute point, and liquid 4 He at the A tran- 
sition point [1]. The boundaries influence the system to a depth given by 
the bulk correlation length £,oo{T, • • • ) (here • • • stays for the dependence on 
other essential parameters, e.g. an external magnetic field H). In the vicinity 
of the bulk critical point (,oo(T) ~ \T — T c \~ u , where v is its critical exponent. 
When £oo(^) <C L the Casimir force, as a fluctuation induced force between 
the plates, is negligible. The force becomes long-ranged when ioo{T) diverges 
which takes place near and below the bulk critical point T c in an 0(n), n > 2 
model system in the absence of an external magnetic field [9]. Therefore in 
the statistical-mechanical systems one can turn on and off the Casimir effect 
merely by changing, e.g., the temperature of the system. 

The effect is of particular experimental interest in studies of wetting of a wall by 
binary liquid mixtures close to their critical end point (see, e.g., [10]). Whereas 
the free energy of magnetic films cannot be measured directly, the free energy 
of liquid films is of relevance for surface tension measurements and wetting 
phenomena. In wetting phenomena thin films are formed such that at least one 
confining boundary is not rigid but determined by thermal equilibrium [7,10]. 
The Casimir force enters into the force balance and thus forms the equilibrium 
thickness of the wetting films, which is accurately measurable. 

The temperature dependence of the Casimir force for two-dimensional systems 
is investigated exactly only on the example of Ising strips [4]. The upper crit- 
ical temperature dependence of the force in 0(n) models has been considered 
in [5]. The only example where the force is investigated exactly as a function 
of both the temperature and the magnetic field scaling variables is that of the 
three-dimensional spherical model under periodic boundary conditions [9]. 
There exact results for the Casimir force between two walls with a finite sep- 
aration in an L x oo 2 mean-spherical model has been derived. The force is 
consistent with an attraction of the plates confining the system . The most 
results available at the moment are for the Casimir amplitudes. For d = 2 by 
using conformal-invariance methods they are exactly known for a large class 
of models [1]. In addition to the flat geometries recently some results about 
the Casimir amplitudes between spherical particles in a critical fluid have been 
derived too [11]. For d ^ 2 results are available via field-theoretical renormal- 
ization group theory in 4 — e dimensions [1,5,11], Migdal-Kadanoff real-space 
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renormalization group methods [12], and, relatively recently, by Monte Carlo 
methods [13]. 

In addition to the statistical mechanics the Casimir forces are object of in- 
vestigations also in the quantum electrodynamics. Nowadays the effect is also 
presented in studies of topics like nonplanar geometries of conducting or dielec- 
tric macroscopic bodies immersed in fluids, the bag model for the description 
of quark confinement in hadrons according to the quantum chromodynamics, 
Casimir effect in general curved space-time in cosmology, models of the early 
Universe, etc. For a review on the Casimir effect in the aforementioned fields 
one can consult, e.g., [2]. A relatively recent review on the Casimir effect in 
statistical-mechanical systems can be found in [1]. 

In recent years there has been a renewed interest [14,15] in the theory of 
zero-temperature quantum phase transitions. Distinctively from temperature 
driven critical phenomena, these phase transitions occur at zero temperature 
as a function of some non-thermal control parameter (or a competition be- 
tween different parameters describing the basic interaction of the system), and 
the relevant fluctuations are of quantum rather than thermal nature. 

It is well known from the theory of critical phenomena that for the temper- 
ature driven phase transitions quantum effects are unimportant near critical 
points with T c > 0. It could be expected, however, that at rather low (as 
compared to characteristic excitations in the system) temperatures, the lead- 
ing T dependence of all observables is specified by the properties of the zero- 
temperature critical points, which take place in quantum systems. The dimen- 
sional crossover rule asserts that the critical singularities of such a quantum 
system at T = with dimensionality d are formally equivalent to those of a 
classical system with dimensionality d+z (z is the dynamical critical exponent) 
and critical temperature T c > 0. This makes it possible to investigate low- 
temperature effects (considering an effective system with d infinite space and 
z finite time dimensions) in the framework of the theory of finite-size scaling 
(FSS). The idea of this theory has been applied to explore the low-temperature 
regime in quantum systems [14-16], when the properties of the thermodynamic 
observables in the finite-temperature quantum critical region have been the 
main focus of interest. The most famous model for discussing these properties 
is the quantum nonlinear 0(n) sigma model (QNLcrM) [14,16]. 

Let us note that an increasing interest related with the spherical approximation 
(or large n-limit) generating tractable models in quantum critical phenomena 
has been observed in the last few years [17-22]. There are different possi- 
ble ways of quantization of the spherical constraint. In general they lead to 
different universality classes at the quantum critical point [17-20]. 

In this paper a theory of the scaling properties of the free energy and the 
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Casimir forces of a quantum spherical model [17] with nearest -neighbour and 
some special cases of long-range interactions (decreasing at long distances r 
as \/r d+a ) is presented. Only the film geometry L x oo d_1 x L T (where L T ~ 
h/iksT) is the finite-size in the imaginary time direction) will be considered. 
The plan of the paper is as follows: we start with a brief review of the model 
and the basic equations for the free energy and the quantum spherical field in 
the case of periodic boundary conditions (Section 2). Since we would like to 
exploit the ideas of the FSS theory, the bulk system in the low-temperature 
region is considered like an effective (d+z) dimensional classical system with z 
finite (temporal) dimensions. This is done to make possible a comparison with 
other results based on the spherical type approximation, e.g., in the framework 
of the spherical model and the QNLaM in the limit n — > oo. The scaling forms 
for the excess free energy, the spherical field equation and the Casimir force 
are derived for a |<r < d < |<r dimensional system with a film geometry in 
Section 3. In Section 4 we present some results for the Casimir amplitudes in 
the case of short-range interactions and in some special cases of long-range 
interactions. The paper closes with concluding remarks given in Section 5. 



2 The model 

The model we will consider here describes a magnetic ordering due to the 
interaction of quantum spins. It is characterized by the Hamiltonian [17] 



where Se are spin operators at site I. The operators Vt play the role of "con- 
jugated" momenta (i.e. [S e ,S e ,] = 0, [V £ ,V e ,] = 0, and [V e ,S e ,] = iS u ,, with 
h = 1). The coupling constant g measures the strength of the quantum fluctu- 
ations (below it will be called quantum parameter), H is an ordering magnetic 
field, and the spherical field /i is introduced so as to ensure the constraint 



Here N is the total number of the quantum spins located at sites of a finite 
hypercubical lattice A of size L\ x L2 x • • • x Ld = N and (■ • • ) denotes the 
standard thermodynamic average taken with the Hamiltonian Ti. In (4) the 
coupling constants J u , are decreasing at large distances — as l/\£ — £'\ d+a , 
where a determines the range of the interaction: i) < o < 2 for long-range 
interaction and ii) a > 2 for short-range interaction. 

The free energy of the model in a finite region A under periodic boundary 



(4) 



ft = H - \ J w s t s t> + \v H s l - H 5 < 




(5) 
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conditions applied across the finite dimensions has the form [22] 



(6) 



PU ((3, g, tf) = sup{i£ In [2 sinh (±/3u; (q; /.))] - |/3 - J^^ 2 } . 

Here the vector qr is a collective symbol, which for Lj odd integers has the 

components j^ 1 , • • • , ^rf \ > n j e { — "^r^ ' ' ' ' ^F^}' an< ^ ^ * s ^ e i nverse 
temperature with the Boltzman constant k B = 1. In (6) the spectrum is 
(q; /i) = g (// + f/(qr)) with £/(<z) = 2JJ2i=i (1 — cosgj) for nearest neigh- 
bour interactions and will be taken of the form i7(qr) = J p a \q\ a , < a < 2, 
for long range interactions (p a > is a parameter to be taken equal to one in 
the remainder). In the above expressions U(q) is the Fourier transform of the 
interaction matrix where the energy scale has been fixed so that U(0) = 0. 
The supremum in Eq. (6) is attained at the solutions of the mean-spherical 
constraint, Eq. (5), that reads 

t 00 1 h 2 



iv m ^ oo Y^ + f/ (9)/ J + ^ 



2 ' A2- 



where we have introduced the notations: b = (2irt) /A, A = y g/J is the normal- 
ized quantum parameter, t — T/J- the normalized temperature, h = H/y/j 
- the normalized magnetic field, and <fi = /i/J is the scaled spherical field. 
Eqs. (6) and (7) provide the basis of the study of the critical behaviour of the 
model under consideration. 



In the thermodynamic limit it has been shown [17] that for d > a the long- 
range order exists at finite temperatures up to a given critical temperature 
t c (\). Here we shall consider the low-temperature region for \a < d < |<r. 
We remind that |<j and |<r are the lower and the upper critical dimensions, 
respectively, for the quantum critical point of the considered system. 



3 Scaling form of the excess free energy and the Casimir force at 
low temperatures 



For a system with a film geometry L x oo d_1 x L T (where \a < d < |<r), after 
taking the limits L 2 — > oo, • • • , — > oo in Eq. (6) and by using the Poisson 
summation formula for the only remaining finite space dimensionality L\ = L, 
we receive the following expression for the full free energy density 
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A k d ( d\ d + i ~ / A A / A 



£L 2 ™ r dx d T . T Nl 
(2vrj2 m=1 ^ m 2 



2sinh( — \/o - - ,v n 



where /c^ 1 = |(47r) ar(d/2), is the radius of the sphericalized Brillouin 
zone, K u (x) and J„(rr) are the MacDonald and Bessel functions, respectively, 
and the critical value of A = A c is 

(9) K l = \{2^Y d jd d q {U{q)/J)^. 

In Eq. (8) <fi is the solution of the corresponding spherical field equation that 
follows by requiring the partial derivative of the r.h.s. of Eq. (8) with respect 
to <fi to be zero. The bulk free energy fbuik(t,\,h) results from f(t,\,h;L) 
by merely taking the limit L — > oo in it. Let us denote the solution of the 
corresponding bulk spherical field equation by (poo. Then for the excess free 
energy it is possible to obtain the finite size scaling form 

(10) f ex {t, A, h; L)/L = XL~^X{x 1 , x 2 , a), 

with scaling variables x\ = L^ 1 ^ (1/A — 1/A C ), x 2 = hL A ^ and a = tL z /\. 
Here v~ x = d — |er, Ajv — | (d + |<t) and z — \a are the critical exponents 
of the model [17]. In Eq. (10) the universal scaling function of the excess free 
energy may obtained in an explicit form (this will be presented in a subsequent 
article) . 

For the Casimir forces in the considered system 

(11) Fc as imir(T, A, h] L) = \L~( d+Z ^ X C asimir{xi, x 2, 0>), 

where the universal scaling functions of the Casimir force X casimir ( x i, x 2, a ) is 
related to that one of the excess free energy X = X(x±, x 2 , a) by 

. . 1 dX A dX OX 

(12) Xcasimir{ x i, x 2,a) = -{d + z)X x ± - 1 x 2 - h za— . 

v ox i v ax 2 da 



4 Casimir amplitudes 



In this section we determine the Casimir amplitudes of the model for the case 
of short range interactions at d — 2 and for some special cases of long range 
interactions. 
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4-1 short range interactions (d = 2) 



It can be shown that the solution y of the spherical field equation for the 
finite system with a film geometry L x oo x L T at zero temperature (i.e. 
1 <C L < oo, L T = oo) is yo = In (\/5/2 + 1/2^ at the quantum critical 
point A = A c , h = (at this point = 0) [21,22]. Setting this value of yo 
in the scaling function (10) of the excess free energy, taking into account that 

1/x), the identity In 



Kz 

2 



7r/ (2x) exp(— x)(l 



2sinh ( y/y/2 



Vy/2- 



Lii expy—y/y) and the properties of the polylogarithm functions Li p (a;) [23], 
we obtain from Eqs. (3) and (10) that the Casimir amplitude is 



(13) 



Aperiodic b.c. 



203) 

57T 



-0.153051. 



Here ((3) is the Riemann zeta function. 



4-2 long range interactions (d = a) 



Let us consider the "temporal Casimir amplitude" in a system with a geometry 
oo d x L T with d/a — 1 at the quantum critical point A = A c , h = 0. In a way, 
similar to that one explained for the case considered previously one obtains 
(0 < a < 2) 

16 t(3) t 3 
(14) /(f,A c ,0;oo)-/(0,A c ,0;oo) = -A c ' ^" > 



: 5a (47r)-/ 2 r( ( r/2)' 
From here one can identify the "temporal Casimir amplitude" to be 

16 C(3) 1 



(15) A t (a) 



5a (47r) CT / 2 r( ( j/2)' 



5 Concluding remarks 



In the present article the free energy of a system with a geometry L x oo^ 1 x L T 
(where \a < d < |<r), is derived (see Eq. (8)). For a = 2 this new result 
reduces to the one reported in [22] where only the case of short-range inter- 
actions has been considered. A general expression for the Casimir force in the 
quantum spherical model is obtained (see Eqs. (11,12). In the classical limit 
(A = 0) for a system with short-range interaction it coincides with the corre- 
sponding one derived in [9] for the classical spherical model. Except for some 
two-dimensional systems the only model within which the Casimir effect is 
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investigated in an exact manner for d 7^ 2 is that one of the spherical model. 
The investigations presented here complement the aforementioned ones for the 
case when quantum fluctuations are present in the system. In order to derive in 
a simple closed form the Casimir amplitudes some particular cases have been 
considered (d — a). For d — 2 this amplitude is given in Eq. (13). This am- 
plitude is equal to the "temporal Casimir amplitude" (i.e. the corresponding 
temperature corrections in an oo 2 x L T system to the ground state of the bulk 
system) for the 0(n) sigma model in the limit n — > 00 [23]. We have demon- 
strated here explicitly that the two models, due to the fact that they belong to 
the same universality class, indeed possess equal Casimir amplitudes as it is 
to be expected on the basis of the hypothesis of the universality. We note that 
instead of considering a finite system with a film geometry at the zero tem- 
perature one can consider a bulk system at low temperatures. As it is already 
clear from above this leads to the same result f or the Casimir amplitudes. 
Note, that in accordance with the general expectations these amplitudes are 
negative. The correction to the ground state energy of the bulk system due to 
the nonzero temperature is determined for the general case d = a by Eq. (15). 
Note, that the defined there "temporal Casimir amplitude" A t (a) reduces for 
a = 2 to the "normal" Casimir amplitude, given by Eq. (13). This reflects the 
existence of a special symmetry for that case between the "temporal" and the 
space dimensionalities of the system. 
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